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1 Introduction 

In a recent paper (cf. |3|) the authors have investigated 3-body motions with 
vanishing angular momentum, in the framework of equivariant Riemannian ge- 
ometry and by resuming the basic approach dating back to Jacobi's geometriza- 
tion of Lagrange's least action principle, in the setting of kinematic geometry 
of 3-body systems. A geometric reduction method was described which reduces 
the study of trajectories of 3-body motions, first from the level of the configura- 
tion space to the level of the moduli space of congruence classes of m-triangles, 
and then 
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• further reduces the moduh curves to that of their shape curves on the 
2-sphere S"^. Namely, a trajectory of 3-body motions is completely deter- 
mined, up to global congruence, by its shape curve which only records the 
changing of shape (i.e. similarity class). 

• Moreover, the unique parametrization theorem further proves that the 
trajectory is already determined by the geometric (i.e. non-paramctrized) 
shape curve. 

• Another remarkable property of the above shape curves is expressed by the 
monotonicity theorem, concerning their (mass modified) latitude function 
on the sphere S^. 

The monotonicity theorem is definitely only valid in the case of zero angular 
momentum, but with this paper we shall extend the first two of the above three 
stated results to the more general case of planary motions. We start with a 
description of the basic setting from |3] in the two subsections below, and a 
summary of the major results is presented in Section 1.3. 

In Section 2 we work out the kinematic and dynamical metric on the moduli 
space, together with the associated differential equations. Finally, in Section 
3 we establish the remaining results needed for the proofs of the two main 
theorems stated in Section 1.3. 



1.1 The basic kinematic quantities and the potential func- 
tion 

The classical 3-body problem in celestial mechanics studies the local and global 
geometry of the trajectories of a 3-body system, namely the motion of three 
point masses (bodies) of mass > 0, say normalized to '^rui = 1, under the 
influence of the mutual gravitational forces. This system constitutes a conser- 
vative mechanical system with the Newton's potential function 

-.. = |a.-a,| (1) 

and potential energy —U. We introduce the vector b = (ai,a2,a3), called an 
m-triangle, which records the position of the system in an inertial frame with 
the origin at the center of mass, and hence ^ m^ai = 0. 

A trajectory is a time parametrized curve ^{t) representing a motion of the 
3-body system, locally characterized by Newton's equation 

, , I dU I dU I dU , 
7 = VC/(7) = (— — — ) 2 

dt^ mi oai m2 c'a2 oa^ 

However, the trajectories can also be characterized globally as solutions of a 
suitable boundary value problem, characterized as extremals of an appropriate 
least action principle, such as the two principles due to Lagrange and Hamilton. 
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Let us also recall the basic kinematic quantities which are the (polar) moment 
of inertia, kinetic energy and angular momentum, respectively defined by 

/ = ^mi|aj|^, T^^^m,\aif, ft ^ ^mi{a, x a,) (3) 

The dynamics of the 3-body problem is largely expressed by their interactions 
with the potential function U via the equation and for example, the invari- 
ance of the total energy 

h = T-U (4) 

is a simple consequence of (jSJ and the definition of T. On the other hand, 
whereas the invariance of the vector O follows from the rotational symmetry of 
U, in this article we shall exploit the consequences of the additional homogeneity 
property of U. 

1.2 Reduction to the moduli space and the shape space 

In this article we shall only be concerned with planary three-body motions, 
namely the m-triangles 6 are confined to a fixed plane and hence belong to 
the configuration space 

3 

1=1 

With the inner product of m-triangles defined by 

5 ■ S' = m^ai ■ hi (6) 

M is given the kinematic metric^ namely the metric such that the right side 
of Newton's equation Q is the gradient of U . Then the squared norm is the 
moment of inertia, / = I{6) — \5\ , and hence the hyperradius p ~ \fl is the 
natural scaling function which also measures the distance from the origin. 

The rotation group S0{2) acts naturally, and by orthogonal transformations, 
on M by rotating m-triangles, and the orbit spaces of M (resp. its unit sphere 
M-^ ~ S^) are the (congruence) moduli space M (resp. the shape space M*), 
namely 

M = M/SO{2), M* = M^/S0{2) (7) 

The points in M represent congruence classes S of m-triangles, and points in 
M* represent the shapes (or similarity classes) S* of non-zero m-triangles. 

Geometrically,the above orbit space construction and orbit map M M 
is, in fact, just the Hopf map construction, whose restriction — > S'^ is the 
classical Hopf fibration. Namely, the spaces fit into the following diagram 

M ~ R4 — > M 
U U U (8) 

^M*~S'2(f/2) 
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where M ~ R"* is a chosen 5'0(2)-equivariant isometry (i.e. choice of Jacobi 
vectors). In particular, M « K'^ is a cone over M* and there is the radial 
projection M — {O} M* which "reduces" a non-zero congruence class 5 to its 
shape 5* . Note, however, the representation of the various shapes of m-triangles 
on a fixed model sphere S"^ depends on the mass distribution {mi}, via the mass 
dependence of the Jacobi vectors. 

Briefly, in this article we shall analyze the two-step reduction 

Af^M, M-{0} ^Ar, lit) ^ lit) -^i*it) (9) 

by which a trajectory lit) of a planary 3-body motion is projected to its moduli 
curve lit) and further to its shape curve 7*(i) on a 2-sphere. In Section 2.1 
we shall put the above reduction and the spaces involved in the framework of 
Riemannian geometry, and moreover, explain how Jacobi's geometrization idea 
can be reduced and extended to the level of M. 

1.3 A summary of the main results 

The Hopf map construction JSJ makes it convenient to use a Euclidean model, 
M — M."^, for the moduli space and with the unit sphere 5'^(1) as the shape 
space Af*. In this way one can express all kinematic quantities and dynamical 
equations in terms of spherical geometry and spherical coordinates, and hence 
take the full advantage of the cone structure of M over M*. 

One can start from Newton's equation ^ for planary m-triangles 

^7 = Vf/(7), lit)€M (10) 

with any 50(2)-invariant potential function U, and hence it is a function on M. 
The additional crucial property of U that we have exploited is its homogeneity, 
namely it is of type 

U = (11) 

where U* denotes the restriction of U to the sphere M*. The Newtonian case 
e = 1 is certainly the most important one, but the proofs are essentially the 
same for other (integral) values of e > 0. 

Consider trajectories lit) of (|10|l for a given energy-momentum level ih,uj), 
Lu = and for spherical coordinates (p, ip, 9) in M, let the curves 

iit) = ipit),i*it)), i*it) = i^it),eit)) (12) 

be the associated moduli and shape curve, respectively. In Section 2.2 we show 
the reduced Newton's equation in M can be presented as the pair 

i =2iU + 2h) + 2il-e)U, 7* + P7* + gVC/* = (13) 

where the first equation in H13() is simply the Lagrange- Jacobi equation, and 7* 
is the covariant acceleration of 7* as a spherical curve. Moreover, the energy 
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integral |0J is the following first order equation in M 



where v = I'y*] — y ip + (sin ip)6 is the speed of the shape curve. In fact, 
combined with 114|) any of the three scalar equations in can be derived 
from the other ones. The equations of H13|l are presented in the coordinates 
(p, (p, 6) in Section 2.2. 

On the other hand, let K* be the geodesic curvature of 7* and U* the 
directional derivative of U* normal to 7*. Then there is the formula 

p2+e = ^e, where & =^ (cf. ISH) (15) 

which separates the radial variable p from the spherical ones. Using H15|) . the 
dependence on p in the coefficient functions P, Q in (|13() can be eliminated, 
namely 

p~2 + ee 2 + ev-' p2+e " g' l-^'^J 

which yields a third order equation for 7* which is, in fact, independent of (/i, uj). 

Clearly, the above function 6 depends only on the relative geometry between 
7* and the gradient of U* . However, we regard (3 as undefined if 7* is a 
geodesic arc (and hence lies on a gradient line), and any such solution of ifT^ 
is called exceptional. Now, assuming (for simplicity) that the shape curve is not 
of exceptional type, our main results can be formulated neatly as the following 
two theorems : 

Theorem 1.1 For a given total energy and nonzero angular momentum, a pla- 
nary three-body motion is completely determined up to congruence by its time 
parametrized shape curve (which only records the changing of shape). 

Theorem 1.2 The time parametrization is uniguely determined by the relative 
geometry between the oriented geometric (i.e. non-parametrized) shape curve 
and the gradient vector field of U* . 

Remark 1.3 (i) In the case of non-zero total energy, Theorem remains 
unchanged in the case of zero angular momentum, whereas the motion is deter- 
mined up to congruence and scaling in the case of (h,u}) = (0,0). We refer to 
l^. Section 4.2. 

(a) Unigueness of time parametrization means, of course, modulo time trans- 
lation, or modulo an affine time transformation when {h^oj) = (0,0). 

(Hi) The proofs of the above theorems are the same for any homogeneous 
potential function of type with e > 0, and e integral in Theorem M.'A We 

choose the most important case, e = 1, in Section 3.2 and work out the crucial 
details. However, the formulae are even simpler in the case e — 2. 
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2 Riemannian geometry and reduction to the 
moduli space 

2.1 Riemannian structures on the moduli space M 

In his famous lectures ^ , Jacobi introduced the concept of a kinematic metric 
ds^ on the configuration space M of a mechanical system with kinetic energy 
T. For example, in the case of an n-body system with total mass '^mi = 1, 

ds^ = 2Tdt^ = m,{dxj + dyf + dzj) (17) 

i 

which is clearly equivalent to the definition ©. Now, for a system with potential 
energy —U and a fixed total energy ft,, set 

Mh = {p e M; h + U{p) > 0} (18) 
dsl = [h + U)ds^ 

where ds\ is called the dynamical metric on M^. By writing 
dsh = VhTUds = VTds = V2Tdt 

Jacobi transformed Lagrange's action integral (on the left side of (|19|l 'l into an 
arc-length integral, namely 

J{^)^ Jjdt^^l^dsh (19) 

and hence the least action principle becomes the following simple geometric 
statement : 

" Trajectories with total energy h are exactly those geodesic curves (20) 
in the space M/j with the dynamical metric ds^ " 

Nowadays, the metric spaces {M, ds'^), {Mh, dsf^) are called Riemannian 
manifolds, and the dynamical metric is a conformal modification of the kine- 
matic metric by the scaling function ([/ + h). As exemplified by (|17|l . a Rieman- 
nian metric on a manifold N amounts to the choice of a kinetic energy function 
on the tangent bundle, T : TN —^ R, which is a positive definite quadratic form 
on each tangent plane TpN. This allows us to define the arc-length function 
u{t) and the kinetic energy along a given time parametrized curve T(t) by 



1 ,du^^ 1 



dr 
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(21) 



Now, let us determine the appropriate kinetic energy T on the moduli space 
M and hence also define its kinematic metric using the recipe (|21|l . At the 
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same time, referring to the diagram © and the Hopf map, let us also introduce 
the orbital distance metric ds^ on M as an S'0(2)-orbit space. Then {M,ds^) 
inherits the structure of a Riemannian cone over the shape space (M*,dcr^), 
namely 

ds^ = dp^ + p^da^, da^ = ds^lw (22) 

Moreover, it is well known that the Hopf fibration in the above Riemannian 
setting is 5'^(1) S'^(l/2), and consequently 

(M*,dCT2) ~5'2(l/2) (23) 

is also the round sphere of radius 1/2. 

Consider a curve ^{t) in M and its orthogonal velocity decomposition 7 = 
A^h _|_ ^uj corresponding splitting of kinetric energy 

T = t'' + T'^, (24) 

where 7" is tangential to the S'0(2)-orbit, and hence T'^ is the kinetic energy 
due to purely rotational motion of m-triangles. By definition of the metric 
ds^, the orbit map M — > Af is a Riemannian submersion and hence maps the 
"horizontal" component 7'' isometrically to the velocity vector of 7. This shows 
T'' = T is also the kinetic energy at the level of M, that is, the kinematic metric 
coincides with the orbital distance metric, and by H22|) , (|23|l and H24|l the latter 
can be finally expressed as 

ds^ = 2Tdt^ = 2(T - T'^)dt^ = 2(r :^)dt^ (25) 

2p 

= dp^ + p^da^ = dp^ + -^{dif^ + sin^ tpdO^) 

Remark 2.1 The expression dip^ + siv? (pdO^ in the last line of \25]) is the 
metric of the unit sphere S'^(l) in terms of spherical polar coordinates. In fact, 
the metric ds^ differs from the Euclidean metric only by the factor 1/4 in f25\}. 
which makes it singular at the origin. Moreover, ds^ is actually a conformal 
modification of the Euclidean metric (cf. Section 2). 

Next, we turn to the construction of the dynamical metric on the moduli 
space M, which depends on U and a given energy-momentum level (li,uj). Fol- 
lowing the geometrization idea of Jacobi, we want the geodesies of this metric to 
be the trajectories in M, regarded as a simple mechanical system with kinetic 
energy T, potential energy U , and conserved total energy h = T — U . Thus we 
introduce the reduced potential function on M 



and define the dynamical metric 



dsl ^ = fds^ = {U + h)ds^ = {U + h ^)ds2 (26) 

2/9^ 
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Finally, it is not difficult to see that Lagrange's least action principle p9|l as 
well as Hamilton's least action principle using the Lagrange function L = T+U, 
can be pushed down to the level of M. This yields the function L — T + U, and 
for example, by following Jacobi's geometrization idea applied to Lagrange's 
action integral in M 

J(7) = V2 Tdt ^ V2 {tJ + h)dt = / VtJ + hds = / dsh,uj , 
J y J ^ J ^ J 

we arrive at the following geometric statement similar to (|20|l : 

"Curves in M representing trajectories in M at a given 
energy-momentum level (/i, uj) are exactly those geodesic (27) 
curves in M with the induced dynamical metric ds\ j^." 



2.2 The geodesic equations of the moduh space 

The moduli space M is, first of all, equipped with the kinematic metric 




and for each energy-momentum level {h, uj) there is the following conformal 
modification of ds^ 

dsf^_^^ = fds^ ^{U + h)df (29) 

The latter is the dynamical metric which characterizes those moduli curves 
7(t) representing trajectories 7(i) at the specified level {h,u!). Namely, 7 is 
a geodesic of the Riemannian metric (I29II . which in the spherical coordinates 
(p, If, 9) expresses as 

dsl,.) = (^^^ + h ^)[dp' + ^{dp' + sin^ ^d0')] 

The standard procedure for the derivation of the geodesic equations, via the 
calculation of the corresponding Christoffel symbols, yields the following system 
of equations expressed with respect to time as the independent variable : 

(z) = p+^-i(^L/* + 2/i) 
P P P 

(zz) 0= ^ + 2^^-isin(2^)^'-^C/; (30) 

(m) = 9 + 2-9 + 2cot{ip)Lpe ~ ^ ^C/g 

p p +^ sin ip 

Note that equation (i), associated with the radial variable p of Af, as a cone 
over the sphere M* = S*^, is simply the Lagrange- Jacobi equation, cf. 
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Moreover, the second equation in is merely a reformulation of equation (ii) 
and (iii), as explained in Section 3.4.2. On the other hand, the dependence 
on (jj in the above equations is only implicit, but it appears in equation (i) via 
substitution of the energy integral 

h^f-U^\-p'+ P^{(p^ + sin^ ^■e') + (31) 

once we have specified the value of h. Equation (iv) makes any of the three 
equations of (|Sn|l superfluous and may be replaced by (iv), as the first step of 
integration, with w appearing as an integration constant. 

Let us also describe another approach to derive the ODEs in (PUJ, namely 
by regarding M as a simple conservative mechanical system with the Lagrange 
function 

L^f + U = T + U-4=^^p' + ii^' + ^^^'^0')-f, + ^ (32) 

Then, straighforward calculations of the associated Lagranges's equations 

d dL dL d dL dL d dL dL 
Je~d}i'^~d^' di^d^'^d^' Jt^^'^~de ^ ' 

yield the system ()3U|I. Similar calculations are worked out in J], Section 3.2. 



3 The proofs of Theorem 11.11 and Theorem 11.21 
3.1 Separation of the scaling variable 

Let 7*(i) — {ip{t),6{t)) be a given time parametrized curve on the unit sphere 
5^, and set s = s(t) > to be its arc- length function. Then its unit tangent 
and positively oriented unit normal are, respectively 

^ dj* I,, d -9, ^ 1/;,. d ■ 1 

ds V oip Ou V Oifi sm ip at) 

and its speed and scalar acceleration are, respectively 

I ^ d 1 -2 

V = y (p^ + SIT? ipO , V ^ —V = -[ipip + (sm if cos p)Lp9 + sin^(iy9)^^] (35) 

dt V 

One way to calculate the geodesic curvature function K* is to express 7* in 
Euclidean coordinates as x(s) = (x(s), y{s), z{s)) and use the formula 

K*{s)=x{s) X x'(s) •x"(s) 

where differentiation is with respect to s. Then, by returning to spherical coor- 
dinates 

K* = (cos ip)e'{l + ip'^) + sin p{p'e" ~ d'lp") 

= ^^{cos ip)0{v^ + (p'^) + sm(p{(pe - Op)^ (36) 
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Next, let us eliminate the second order terms ip and 9 in the expression 
using equations (ii), (iii) of the system H30|l . This procedure yields 



K*v^ = {cos Lp)9{v^ + <i?) + {smLp)Cp ( -^6 - 2{cot ip)(pe + ^ \—Ug 

\ p p^^^ sin ip 

- {sm^)d (^-^^ + isin(2^)^' + ^C/^ 



p2+e \^si2^(^ p2+e 

and consequently we arrive at the formula 

P'^' = 4^ = 4© (37) 

Note that the function 6 = 6(7*), called the Siegel function in depends 
only on the intrinsic geometry of the pair (7*, U*) on the sphere. 



3.2 Intrinsic geometry of the shape curve and the gradi- 
ent of U* 

In the local analysis of the moduli and the shape curve, and their interaction 
with the potential function U* , we shall distinguish between two types of vari- 
ables or quantities associated with a given moduli curve ^{t) — {p{t) ,j*{t)). 
On the one hand, the intrinsic quantities depend only on 7* as an oriented 
geometric (i.e. unparametrized) curve and U* as a function on S^, and on the 
other hand, the variable quantities are defined along 7 or 7*, depending on the 
scaling function p in the moduli space M or the time parametrization of the 
curves. 

The basic intrinsic quantities are the gradient field VC/* (or its tangential 
and normal derivatives U*, U*), the orthonormal frame field {t*,^} along 7*, 
and the geodesic curvature function K* of 7* . In general, the linkage between 
7* and U* is neatly encoded into the intrinsic function S — U*/K*, introduced 
in H37|l . so we shall assume 7* is not confined to a geodesic circle (in which case 
6 is undefined). 

We choose a (generic) point Pq on 7* , and let s be the arc- length parameter 
of 7* in the positive direction starting from Po-Then the coefficients of the 
following power series expansions 



K* 


= Ko- 


^Kis 


+ K2S^ 


U* 


= Mo -1 


- UlS + 


U2S^ + 


u; 


= Ul A 


- 2U2S ^ 




u* 




- UJlS -\ 


- UJ2s'^ 4 


& 


= 60- 


f Gis - 
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yield intrinsic quantities (or geometric data) localized at the point Po-The co- 
efficients &n are expressible as rational functions of Ki and LUi, and generally, 
let us say the order of a coefficient in (|38|l is the highest order of derivatives of 
local coordinates in its expression. Thus, we say (/3g, f^nd uq are the intrinsic 
geometric data of order at Pq, and for example, ui,u!q and T*\pg have order 
1, and Wn, Un+i (resp. Kn, &n) have order n + 1 (resp. n + 2). 

Let {(fi, 9) denote spherical polar coordinates so that Pq is different from any 
of the "poles" = or tt. We shall expand the coordinate functions of ^{t), as 
well as U* and its partial derivatives, as power series with respect to t : 



p 


= Po H 


^ Pit- 


1-^2*^ + ^3*^ + 




= fo - 


f f^t 


+ ip.2t^ + ip^t^ + ... 


e 


= 00^ 


-Oit- 


F 621^ + 63.1"^ + 


V 


= vo-\ 


- Wl + 


V2t + V2t^ + 


U* 


"0 - 


\- Uit - 


-f U2f + U^t^ + .... 


u; 


= Mo - 




+ fi2t^ + fi^t^ + .... 




= Vo H 







For convenience, some of the initial coefficiens are 

Uq ^U*{ipQ,9o), ui = ^Qipi + t]q9i, etc. 

fo = sin(2(^„), /i = 2cos(2</Jo)v3;^, etc. (40) 
go = sm^iipg), gi = /o(/3i, etc. 

Vl = — [2(^i</j2 + sin((/3o) cos(93o)V'i^i + 2sin^(V3o)^i^2] 

where vi follows from H35|) . and we also write 

sin(2^) = /o + /it + /2t' + .... 
sin^(93) = .90 + git + g2t^ + .... 

We shall regard /Iq; % as intrinsic data, but they depend on the coordinate 
system, of course. 

Below we shall investigate dependence relations among the coefficients in 
such as Pj, ipp 9k and various other coefficients. Some of them are di- 
rectly expressible in terms of the intrinsic data and hence regarded as constants, 
whereas the others are the variables. 

Definition 3.1 The following list of coefficients in the expansions 

PQ,vo;Pi,fi,di;p2,'P2,02 (41) 

will be referred to as the variables of order < 2. The variables of order n are 
p„, (/3„, 9n when n > 0, and pg, vq are the only variable of order zero. 
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Henceforth, assume the above moduH curve j{t) is a solution of the ODE 
system (|^ - (|^ with e — 1. By inserting the power series into the equa- 
tions (i)-(iv) and applying the method of undetermined coefficients, we arrive 
at the following scheme of recursive relations for the variables of increasing order 
0,1,2..: 

^10 : = 2pgp2 + PqPi - - 

^20 : = 2pg(^2 + ^pIpi^i - ^plfoOl - ipo (42) 
E3o:0 = 2gaple2 + 2.goPoPi6'i + plk^iOi - 4.Vo 
E40 -0^ pIpI + '^pU'A + 9oOi) + - 2moPo - 2Vo 
and in general 

:0 = (n + 2)(n+l)p2p^+2 + 

:0 = (n + 2)(n + l)p3^„+2 + (43) 

E^n : = (n + 2)(n + l)5oPo6'«+2 + 

where the remaining terms are of less order since they involve Pi,iPi, Oi for i < 
n + 2. For example, the order of the coefficients u„, /i„,?7„ in H43(l is n. The 
equations E^n for n > will not be needed and hence omitted since they do not 
lead to additional (algebraic independent) relations. 

Now, let us select some independent and recursive relations from the above 
ones, but first we take the basic identity (|37|l and the expression H35(l for the 
speed in the spherical metric, whose zero order terms yield the two identities : 

Eo : pIvI = 46o (44) 

E'o-.vo^ ^'fil+goOl (45) 

We shall use the symbols Ji, J2 etc. to denote various expressions which are of 
intrinsic type. By using the identities £^10 and E40 can be restated as 

Ei:po{pl-2h) + — ^Ji, Ji=2wo-6o (46) 
Po 

E4 ■■ pIp2 - z — -^4, J4 = 7:(-"o + 60) 

2po 2 

Next, the direction i/jq of 7* at the point ((/?q, Oq) is intrinsic; it is also conve- 
niently represented by the unit tangent vector 



1 d d d d 

The coefficients J^, Jg are intrinsic functions, depending on the coordinate 
system, and they are related by the identity 



4 + 9oJe = 1 (47) 
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Therefore, we adjoin to our list (|46|l the two identities 



E2-^i^ J^vo (48) 
E3 : 9i = Jgvo 

Still, we have not used all zero order relations, namely £'20 and £"30, and 
now we state them as 

E5 ■■ pI'P2 + PoPifi ^ J5, Jb = '^PQ + hJg&o (49) 
Ee : Po^2 + PnPi^i — "^67 J& = — ^— — 2— Jio JeSo 

By continuing this way, we obtain for each n > three new relations with 
leading terms as indicated 





= 


2 

PoPn+2 




= 


PoVn+2 


E^n+i 


= 





(50) 



where the triple p„^2j Vn+21 ^n+2 are the variables of highest order n + 2. 

Claim 3.2 It is possible to solve the above recursive relations for the variables 
fl}) completely in terms of the intrinsic local geometric data in the shape space. 



We proceed as follows. At this point, we observe first that there are alto- 
gether 3n-|-8 variables 

PO; ''^O; Pi, iPl, 6'i;P2i V2i(^2\ Pn+2: 'Pn+2' 

involved in Bn+B recursive relations, and the first eight involve only the variables 
up to order 2. However, Eq,E2 and £3 are obviously algebraic dependent due 
to the identity (|47|) , so we shall search for one more independent relation among 
the variables of order < 2. We expect such a relation to involve local intrinsic 
quantities of order (at least) 3, so a natural approach is to differentiate the 
basic identity (|37() involving the function S. Then, evaluation of the resulting 
identity at t = tg yields 



3^ + 2!^ = J,, J,^^ (51) 



Using the expression in (|40|l for vi we can restate the above identity as 

J7 (52) 



3^ + ± 
Po«o 



1 2 

Vi^2 + -^foVi^i +90O1O2 
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By simple calculation and substitution using some of the previous relations Ei 



1 2 

= <Pi(J5 - pIpiVi) + "^1(^/0^9 Po^o) + ^1(50^6 - goplPiOi) 
Po 

= -4— ©0 + vq [J^ J5 + foJ^Jg&o + goJeJe] 
Po 

and by substitution into H52|) . using the identity PqVq — 4&q and the expressions 
for J5 , Jg in (|49|l , this leads to our new identity 

E[:-^^Js, Js^2eo\J^Po + JeVo)-J7=^{2ui~&i) (53) 

where ui is the tangential derivative U* of U* at Pq, cf. H38|l . 
From the system of algebraic equations 

Eo,Eq, Ei,E[,E2, E-i,Ei, 

we can now solve recursively and thus determine the variables 

P0:1'0,Pl,</'l,^l:P2:'/'2:^2, 

successively in terms of the intrinsic data. In fact, this is obvious from the 
structure of the equations, once we have determined pg, wo, Pij namely using the 
three equations Eo,Ei,E[: 

2 

plvl^Aeo, Poipl^2h) + — = J,, -^ = J8 (54) 

Po PqVo 

It follows that Pq is characterized as a positive root of the following polyno- 
mial of order < 2 : 

2hpl ~ (4j|6o - Ji)Po - c^' = (55) 

The case = is discussed in |2|, Section 4.1, and we recall the three cases 
h = 0,h > 0,h < are characterized by the sign of 4J|6o — Ji, namely the 
cases are 



^o = ^(4J| + l)6o, uo < ^(4J| + l)6o, uo > ^{AJi + l)6o (56) 



Now, assume lu > 0. In the case of h — we have clearly 

^° = J, ^1j|6o ' "°>^(4j| + l)6o (57) 
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For h > the solution must be 

Po = ^ ( (4 J| 60 - Ji) + ^J{Ui&o^Jly+8hcJ^'^ (58) 

but as a further characterization we cannot rule out any of the three types of 
constraints l|56|l . 

In the third case, h < 0, the two roots of equation H55|) would be negative if 
4J|6o — Ji > 0, so the inequahty "< " must hold and consequenly 

Po = ^ ( (4^1 So - Ji ) - V^(4J|6o- Ji)2+8/i^2^ (59) 

"0 - ^(4J| + 1)60 > ojVWI > 

The other choice of root in the formula for pg is ruled out by demanding con- 
tinuous dependence on the parameters, e.g. a; ^ should not imply —>■ 0. 

Finally, with the above value for pg, the system (|56|) yields the following 
intrinsic formulae for the two variables wo, Pi, namely 

z;o-2J%, p^^2JsJ— (60) 
V Pi) V ^0 

3.3 Summary and final proofs 

To complete the proofs of Theorem 11.11 and Theorem 11.21 let us start with a 
curve 7(t) — {p{t), f{t), 9{t)) in M which is the moduli curve of a trajectory j{t) 
of a planary 3-body motion. Then the curve j(t) in M is uniquely determined, 
up to a global congruence, by the curve j{t) and the size |r2| of the angular 
momentum vector. We refer to (|2], Theorem B) for the purely kinematic result 
concerning the general lifting of curves j{t) in M to curves in M. On the other 
hand, by the formula 137(1 . the size function p{t) is already determined by the 
shape curve 7*(i) = {(p{t),9{t)) (7* assumed to be non-exceptional), and this 
proves Theorem ll.il 

In Section 3.2 it is demonstrated that the power series expansion of 7*(t) is 
essentially determined by quantities which depend only on the geometric (i.e. 
unparametrized) shape curve 7* and the relative geometry between 7* and the 
gradient vector field VC/*. By "essential" we mean that the same shape curve 
can only be reparametrized in the trivial way, namely by an affine transformation 
of time, in order to remain the time parametrized shape curve of a (planary) 
3-body motion. In view of Remark 1.3 this completes the proof of Theorem 
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